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ABSTRACT 


The  horse-shoe  vortex  that  forms  around  a  wing-wall  junction  is  experimentally 
investigated.  Laser-Doppler  velocimeter  measurements  of  mean  velocity  and  higher  order 
statistics  including  the  third  order  products  are  reported  for  15  measurement  stations  in  one  plane 
to  the  side  of  the  junction  of  a  3:2  elliptical  nose,  NACA  0020  tail  wing  and  a  wall.  The  approach 
Reynolds  number  of  the  air  flow  based  on  momentum  thickness  is  =5940.  The  outer  layer  vortical 
structure  on  the  down-wash  or  wing  side  results  in  a  thin  boundary  layer  and  lower  turbulence 
intensity  due  to  the  redirected  free-stream.  Lateral  pressure-gradients  cause  separation  on  the 
uplifting  side  of  the  vortex  on  the  wall.  Bimodal  histograms  of  the  w  fluctuating  velocity  occur 
under  the  vortex  core  near  the  wall.  A  vortical  structure  with  higher  vorticity  concentration  and 
opposite  sense  to  the  large  vortical  structure  is  located  below  the  large  vortical  structure.  High 
normal  stress  values  are  obtained  at  the  wing-wall  junction.  In  this  wing-wall  junction  region  a 
vortical  structure  of  the  same  rotation  sense  of  the  large  outer  layer  vortex  forms.  Triple  products 
describe  the  diffusion  processes  for  this  type  of  a  flow.  In  such  a  flow  the  shear-stress  angle 
(SSA)  highly  lags  the  flow-gradient  angle  (FGA)  and  the  turbulent  structure  highly  lags  the  mean 
flow.  The  turbulence  diffusion  is  highly  altered  due  to  the  presence  of  the  large  outer  layer  vortical 
structure. 


TABLE  OF  CONTENTS 


ABSTRACT 


Nomenclature . 

List  of  Tables . 

CHAPTER  I  -  INTRODUCTION 


CHAPTER  II  -  RESULTS  AND  DISCUSSION  .  3 

11-1.  Mean  Velocities .  3 

II-2.  Normal  Stresses .  4 

II-2.  Shear  Stresses  .  4 

CHAPTER  III  -  PROBABILITY  DENSITY  FUNCTIONS  OF  u,  v,  w  FLUCTUATIONS  AT 

THE  STATIONS  SURROUNDING  THE  VORTEX  CORE .  6 

III-1.  s/t=-0.247  .  6 

III-2.  s/t=-0.291 .  6 

1 1 1-3.  s/t=-0.336  .  7 

III-4.  s/t=-0.368  . 7 


CHAPTER  IV  •  TRIPLE  PRODUCTS .  8 


CHAPTER  V  ■  DERIVED  QUANTITIES .  12 

V-1 .  Flow  Angle  .  12 

V-2.  Flow  Gradient  Angle  . 12 

V-3.  Shear  Stress  Angle .  12 

V-4.  {dW/dy)/{t/UJ  .  12 

y-S.TKE/Uj .  13 

V-6.  V^JU^ .  13 

V-7.  Uv/IF .  13 

V-8.  Skewness  of  u .  13 

V-9.  Skewness  of  v .  14 

V-1 0.  Skewness  of  w  .  14 

V-11.  Shear  Stress  Magnitude .  14 

V-1 2.  Townsend’s  Structural  Parameter  A, .  15 


V-1 3.  iLIlL  Structural  Parameter 


15 


V-1 4.  Rotta’s  T  Parameter 


15 


t 


CHAPTER  VI  -  VARIATION  IN  FLOW  PHYSICS  AT 


•  SELECTED  LOCATIONS .  16 

VI-1.  Separation  Location  .  16 

VI-2.  Stations  Surrounding  the  Vortex  Core  .  17 

VI-3.  Near  the  Wing/Wall  Junction .  18 

•  CHAPTER  VII  -  CONCLUSIONS .  21 

REFERENCES .  23 

TABLES .  25 


APPENDIX  A:  Line  Plots  of  Flow  Variables  Products  and  some  Structural  Parameters  ...  51 


% 


III 


NOMENCLATURE 


static  pressure  coefficient 

• 

Oc,/ax)^, 

static  pressure  coefficient  gradient 

Cij 

dujj-  dup 

'  Reynolds’  stress  convection  tensor 

df 

FA 

arctan{W/U),  flow  angle 

• 

FGA 

arctan(  flow  gradient  angle 

aUlay 

P 

fluctuating  pressure  component 

• 

P|jk1 

production 

triple  product  transport 

tensor 

f^ijk2 

production 

-{u.,UjU^  Uj^+UjU,u„  UiJ,  triple  product  transport 

tensor 

PRii 

_ dil:  _ dU, 

-(uiui^+uju,^),  Reynolds’  stress  production  tensor 

uXj  uXj 

PD, 

1  d(pUi)  dCpu,) 

^  ^  Reynolds’  stress  pressure  diffusion  tensor 

p  dX,  dXj 

# 

Re^ 

S.J 

Reynolds  number  based  on  momentum  thickness 

rate-of-strain  tensor 

SSA 

arctan(^),  shear  stress  angle 
uv 

TD, 

d{u,u,u) 

‘  ,  Reynolds’  stress  turbulent-diffusion  tensor 

TKE 

(iF-i^+W)/2  =(  Ui  Ui)l2  ,  Turbulent  kinetic  energy 

Ue 

velocity  magnitude  at  the  boundary  layer  edge 

• 

Wi 

mean  velocity  components 

Reynolds  stress  tensor.  i=1 ,  2,  3;  j=1 ,  2,  3. 

U^UjU, 

triple  velocity  correlation  tensor 

# 

U* 

nondimensional  mean  velocity 

t/x 

sjijp  ,  skin  friction  velocity 

IV 


•  VD, 


V — Li,  Reynolds’  stress  viscous  diffusion  tensor 
dxf 


^TC’yrc  <^TC 
y*-yuj\ 

•  Pfs 

Pivc 

h 

8 

# 


#  e 

V 

p 


♦  % 

«/y 

dijjk 

^iik 

interaction 

^ijk 

#  dissipation  tensor 


tunnel  coordinates 
wall-law  variable 

flow  angle  at  the  boundary  layer  edge 
wall-stress  direction 
Dirac  delta  function.  5,^=  1  if  i=j,  5,^=0  if 
boundary  layer  thickness 


du,  du, 

-2v — 1 — 1 ,  Reynolds’  stress  viscous  dissipation  tensor 

dXj  dXf 


e„/2,  dissipation  of  turbulent  kinetic  energy 


1  -—)—dy,  momentum  thickness 

uj  U. 


kinemaWc  viscosity 
density 

standard  deviation 

^(-pt7v)^+(-pvw)^ ,  shear  stress  magnitude  in  the  flow 
wall  shear  stress 


n  dUi  du, 

i-{ — 1+ — i),  Reynolds’  stress  pressure-strain  tensor 
p  dXj  dx, 

vorticity  tensor 

+(p/p)u,u;5*„ 

triple  product  transport  diffusion  tensor 

(p/p)L( u,Uy)_*+( UyU*)  ,j ,  triple  product  transport  pressure 
tensor 

2  V  ( + UjU^„  u^J ,  triple  product  transport 


V 


Table  1.  Laser-Doppler  Velocimeter  Locations  and  Flow  Parameters . 

Table  2,  21:1  odds  ±2o  Uncertainties  in  Measured  Quantities  . 

Table  3.  Some  Length  Scales  Obtained  From  LDV  Data  in  22.4°  Line  Coordinates 


CHAPTER  I.  INTRODUCTION 


Wing-body  junction  flows  have  been  the  subject  of  many  studies  in  the  past.  The  flow  at 
the  nose  of  the  wing  has  aperiodic  (bimodal)  double-peaked  fluctuating  velocity  histograms  with 
very  large  Reynolds  stress  values  compared  to  the  approaching  boundary  layer  values 
(Devenport  and  Simpson,  1990a).  Shinpaugh  and  Simpson  (1995)  used  a  scanning  LDV  and  a 
pressure  transducer  to  further  investigate  this  velocity  and  the  pressure  field  simultaneously. 

Martinuzzi  and  Tropea  (1993)  studied  flow  around  several  rectangular  shaped  obstacles 
to  show  the  flow  phenomena  around  the  obstacles  in  a  fully  developed  air  channel.  The  approach 
boundary  layer  separates  due  to  the  adverse  pressure  gradient  in  this  region  and  rolls  in  towards 
the  floor,  generating  backflow.  The  flow  visualization  water  tunnel  studies  of  Khan  et  al.  (1 995) 
and  Kim  (1991)  also  show  these  characteristics  and  describe  the  time  development  of  the  motion. 
The  vortical  flow  is  stretched  around  the  wing  and  is  reoriented  by  the  wing,  which  results  in 
vortical  flow  along  the  sides  of  the  wing-body  junction. 

The  3-D  TBL  of  the  pressure-driven  type  has  been  the  subject  of  many  papers.  Reviews 
of  recent  experimental  data  were  done  by  Anderson  and  Eaton  (1989),  Flack  and  Johnston 
(1994),  Schwarz  and  Bradshaw  (1992),  and  OIgmen  and  Simpson  (1992, 1993).  DNS  studies  of 
Spalart  (1989)  and  Moin  et  al.  (1992)  also  show  the  capability  of  CFD  in  explaining  the  flow 
physics  at  low  Reynolds  numbers.  Three-dimensional  boundary  layer/vortex  interaction,  which  is 
relevant  to  the  work  undertaken  here,  was  studied  by  Eaton’s  group  (Eaton,  1995). 

In  the  present  study  the  detailed  mean  flow,  Reynolds  stress  and  triple  product 
characteristics  of  a  wing/body  junction  flow  in  one  plane  at  the  side  of  a  NACA  0020  tail;  3:2 
elliptical  nose  wing  (maximum  thickness  t=  7.17  cm)  are  presented.  Figure  1  shows  the 
measurement  locations  and  the  coordinate  axes  used  in  the  present  study.  All  the  data,  unless 
othenwise  specified,  are  presented  in  a  coordinate  system  where  X22.4  is  perpendicular  to  the  plane 
of  measurements  "s"  is  in  the  plane  of  the  measurements  with  positive  axis  towards  the  wing  and 
y  axis  is  perpendicular  to  the  wall.  Table  1  gives  the  flow  characteristics  at  the  measurement 
locations. 

This  flow  is  extremely  well  documented  and  detailed  references  can  be  found  in  the  work 
by  Simpson  (1995, 1996)  and  OIgmen  and  Simpson  (1995a,  1996e).  Me  Mahon,  Merati  and  Yoo 
(1987),  Me  Mahon,  Hubbart  and  Kubendran  (1982),  and  Dickinson  (1986)  made  measurements 
with  a  similar  model  at  other  Reynolds  numbers  to  report  the  flow  field  around  and  in  the  wake 
of  the  wing  using  hot-wire  probes.  The  work  undertaken  here  is  similar  to  the  work  of  Devenport 
and  Simpson  (1992)  and  of  Fleming  et  al.  (1993)  and  of  the  previously  mentioned  researchers. 
Fleming  et  al.  (1993)  made  measurements  at  several  planes  along  the  side  and  in  the  wake  of 
the  wing  and  described  the  global  flow  characteristics  of  the  junction  vortex  in  comparison  to  the 
data  of  McMahon  (1987)  and  Devenport  and  Simpson  (1990b).  The  work  of  Devenport  and 
Simpson  (1992)  gives  the  general  characteristics  of  the  mean  flow  and  Reynolds’  stresses  at  five 
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different  planes  around  the  nose  and  side  of  the  wing  using  LDV.  They  further  used  their  data  to 
evaluate  some  turbulence  models.  In  contrast  to  the  body  of  existing  experimental  data  this  flow 
has  not  yet  been  extensively  studied  computationally.  Sung  and  Yang  (1988)  showed  that  simple 
algebraic  eddy  viscosity  models  did  not  perform  well  for  such  a  flow. 

The  present  study  differs  from  these  previous  studies  with  more  detailed  and  more 
accurate  near  wall  information  of  not  only  of  the  mean  flow  and  Reynolds  stresses  but  also  of  the 
triple  products.  The  approach  flow  conditions,  static  pressure  distribution  on  the  wall  were 
previously  presented  (OIgmen  and  Simpson,  1990,  1995a).  The  five-component  laser-Doppler 
velocimeter  used  for  the  measurements  and  the  data  reduction  techniques  were  previously 
discussed  (OIgmen  and  Simpson,  1995a,  b).  The  effective  probe  volume  had  a  diameter  of  30 
microns.  The  probe  volume  could  traverse  as  close  as  50  microns  to  the  wall  and  the  wall  could 
be  found  with  an  uncertainty  of  ±10  microns.  The  data  acquisition  and  reduction  are  described 
here  very  briefly.  The  coincident  three  velocity  components  in  the  optical  coordinate  system  were 
transformed  to  the  tunnel  coordinates.  The  skirts  of  the  histograms  of  the  velocity  were 
interactively  cleaned  of  the  noise  using  a  parabola  fitting  routine  to  the  logarithm  of  the 
histograms.  The  inverse-velocity-magnitude  bias-corrected  histograms  calculated  using  the 
cleaned  data  were  used  to  extract  the  mean  velocity  and  higher  order  terms.  Table  2  gives  21 :1 
uncertainties  calculated  using  two  separate  data  sets  measured  at  one  of  the  measurement 
locations  and  using  Chauvenet’s  criterion  to  define  the  uncertainty. 

In  the  following  sections  selected  data  are  presented  and  the  differences  in  flow 
characteristics  at  selected  locations  due  to  different  flow  conditions  are  explained.  Remaining 
quantities  can  be  found  in  the  report  by  OIgmen  and  Simpson  (1996d).  That  report  includes  the 
mean  velocity,  Reynolds’  stress,  all  triple  product  contour  plots,  discussion  of  the  bimodal  nature 
of  the  fluctuating  velocity  probability  density  distributions  and  contour  plots  of  some  derived 
quantities  such  as:  skewnesses  for  different  velocity  components,  turbulent  kinetic  energy, 

longitudinal  vorticity  distribution,  W/IF,  shear  stress  magnitude,  Townsend’s  structural  parameter 

a„  (iF+vP)/i^,  and  velocity  magnitude.  The  transport  equation  balance  for  the  selected  locations 
in  this  flow  were  reported  by  OIgmen  and  Simpson  (1996a,  b),  and  the  evaluation  of  some  of  the 
existing  pressure-strain  and  of  the  turbulent-diffusion  models  for  the  same  stations  of  the  previous 
work  were  reported  by  OIgmen  and  Simpson  (1996a,c  and  1997). 
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CHAPTER  II.  RESULTS  AND  DISCUSSION 


In  this  section  the  mean  velocity  at  15  locations  that  lie  in  a  plane  perpendicular  to  the 
wing  and  the  wall  which  passes  through  station  5  (OIgmen  and  Simpson,  1995a)  are  presented. 
The  X22.4  axis  perpendicular  to  this  plane  lies  at  22.4  degrees  with  respect  to  the  tunnel 
coordinates  and  the  centerline  of  the  wing  (Fig.  1). 

IM.  MEAN  VELOCITIES 

The  approaching  2-D  turbulent  boundary  layer  (TBL)  on  the  wall  separates  from  the  wall 
along  and  around  the  symmetry  line  of  the  wing  due  to  the  adverse  pressure  gradients  generated 
by  the  presence  of  the  wing.  The  separated  shear  layer  further  rolls  in  towards  the  wall-wing 
junction  carrying  the  high  speed  free-stream  fluid  and  generates  a  high  vorticity  flow  zone  which 
is  non-periodic  in  time.  The  vorticity  within  the  approach  2-D  TBL  and  the  vorticity  generated 
within  the  nose  section  of  the  wing  are  stretched  around  the  wing-wall  junction  and  their  direction 
is  reoriented  along  the  wing  contour.  The  corkscrew  motion  of  the  large  rotational  region  shapes 
the  structure  of  the  entire  wing-wall  junction  flow.  The  outer  layer  large  vortical  structure  core  is 
located  at  logio(y/f)*-1.1  and  s/t=-0.32. 

The  high  speed  free-stream  fluid  directed  towards  the  wall  by  the  outer  layer  vortex  results 
in  high  streamwise  U  velocity  components  in  the  zone  between  the  vortex  core  and  the  wing  and 
faster  moving  fluid  as  the  wing  is  approached  (Fig.  2).  (The  boundary  layer  on  the  wing  is  not 
studied  here).  The  dUldz  of  the  flow  is  positive  which  results  in  acceleration  in  the  +s  direction. 

Near  the  wall  dD/dy  also  increases  towards  the  wing.  In  a  region  between  the  vortex  core 
and  the  wing  dUldy  is  approximately  zero  where  a  law-of-the-wall  for  the  U  component  would  not 
hold.  The  boundary  layer  on  the  wall  becomes  thinner  compared  to  the  layer  thickness  at  Station 
5,  which  is  the  station  farthest  away  from  the  wing  presented  in  this  study.  The  free-stream  fluid 
directed  towards  the  wall  also  results  in  high  negative  V  velocity  components  closer  to  the  wing 
and  away  from  the  wall  (Fig.  3),  Positive  V  velocities  are  observed  in  the  uplifting  side  of  the 
vortex.  The  rolling  motion  within  the  measurement  plane  at  locations  closer  to  the  wall  redirects 
the  flow  away  from  the  wing  parallel  to  the  floor. 

This  flow  has  high  speed  negative  W  velocity  components  faster  than  the  outer  layer 
positive  W  free  stream  velocity  component  (Fig.  4).  The  lateral  pressure  gradient  dCpIds  values 
obtained  on  the  wall  in  the  plane  of  the  measurements  are  tabulated  in  Table  1  and  show  adverse 
pressure  gradients  in  the  -s  direction.  The  high  speed  W  velocity  energizes  the  near  wall  flow, 
although  the  adverse  lateral  pressure  gradients  shape  the  lateral  boundary  layer  away  from  the 
wing.  The  near  wall  low  momentum  fluid  affected  by  the  lateral  pressure  gradients  separates  at 
s/t=-0.54  as  shown  by  the  convergence  of  the  stream  lines.  Near  the  wing  and  the  wall 


intersection,  the  W  velocities  change  sign  and  become  positive. 

11-2.  NORMAL  STRESSES 

High  speed  free-stream  flow  directed  towards  the  wall  by  the  large  outer  layer  vortex  also 
carries  the  low  level  turbulence  of  the  free-stream  flow.  This  results  in  low  turbulence  levels 

between  the  vortex  core  and  the  wing  (Figures  5,  6  and  7).  In  this  region  the  the 

dy 

largest  IF  transport  equation  production  term  is  small  since  is  small.  The  IF  maximum  is 

at  a  location  close  to  the  wing  and  wall  intersection  (Fig.  5).  High  u  fluctuations  are  partly  due  to 
redirection  of  the  large  instantaneous  negative  V  velocity  component  to  instantaneous  W  and 
instantaneous  U  velocity  components.  The  v  fluctuations  are  suppressed  due  to  the  presence  of 
the  wall  and  w  fluctuations  are  suppressed  due  to  the  presence  of  the  wing  near  the  intersection. 

High  IF  values  are  also  located  in  the  vicinity  of  the  vortex  core,  between  the  vortex  core  and 

the  wall,  and  near  the  wall  below  the  vortex  core.  Below  the  rotational  separation  region!/^ 
values  are  as  low  as  at  Station  5. 

Maximum  IF  normal  stress  location  is  right  below  the  vortex  core  (Fig.6)  and  highT^ 
values  are  confined  to  a  region  surrounding  the  vortex  core.  The  'F  near  the  wall  is  not  effected 

by  the  outer  layer  vortex.  Both  at  the  wall-wing  junction  and  also  near  the  separation  location ‘F 
values  are  low. 

The  FF  normal  stress  maximum  is  at  a  location  close  to  the  wall  below  the  vortex  core 

(Fig.  7).  The  FF  maximum  is  5/3  times  larger  than  the  IF  maximum.  Large  "F  values  are  due 
to  non  periodic  meandering  of  the  vortex  core  which  will  be  discussed  later  in  this  section.  Similar 

to  the  F  normal  stress  the  FF  normal  stress  is  also  not  affected  by  the  large  scale  mean  flow 

motion  at  the  wing/wall  junction.  The  lowest  "vF  is  located  at  Station  5.  The  IF  values  around 
separation  location  are  close  to  Station  5  values. 

11-3.  SHEAR  STRESSES 

The  velocity  correlations  W,  Uw  and  Tw  are  large  negative  values  both  when  the  two  of 
the  fluctuating  velocities  are  highly  correlated  with  opposite  signs  and  also  when  seldom  occurring 
large  opposite  signed  fluctuating  velocity  couple  outweigh  the  Reynolds’  averaged  mean  value 

towards  large  negative  numbers.  Small  magnitude  Fv  values  are  located  at  the  uplifting  side  of 
mean  vortical  flow  (Fig.  8).  This  region  is  far  away  from  large  F  regions  but  closer  to  the 
maximum  F  location.  At  the  locations  where  the  u  and  v  fluctuations  are  large  either  they  are 
not  of  the  opposite  sign  or  they  are  not  well  correlated.  At  regions  where  Fv  magnitude  is  small, 
the  F  values  are  as  high  as  the  maximum  value  at  Station  5,  however  much  smaller  than  the 
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values  at  the  w^ing-wall  junction  and  than  the  values  below  the  vortex  core.  At  the  same  region 

values  are  twice  as  iarge  of  the  TF.  The  IJv  shear  stress  magnitude  is  much  smaller  at  the 

locations  where  IP  contours  show  maxima.  Positive  Uv  values  are  located  near  the  wing  closer 

to  layer  edge.  In  the  downwash  side  of  the  large  outer  vortical  flow  Pv  shear  stress  magnitudes 
are  smaller  than  that  of  at  Station  5.  High  speed  low  fluctuation  free-stream  fluid  redirected  by  the 
vortex  towards  the  wall  in  this  region  results  in  less  correlated  u  and  v  fluctuating  velocities. 

Around  the  separation  location  the  Pv  stress  values  are  as  low  as  at  Station  5. 

Maximum  value  of  the  Pw  stress  is  obtained  below  the  vortex  core  near  the  uplifting  side 

of  the  vortex  (Fig.  9).  This  location  is  closer  to  the  wall  than  the  Pv  stress  contours  show  a 
minimum.  Another  maximum  occurs  very  close  to  the  wing/wall  intersection.  Minimum  values  of 

the  uw  stress  are  obtained  in  the  downwash  side  of  the  outer  vortical  flow.  These  locations 
roughly  correspond  to  the  locations  where  TP  normal  stress  is  large. 

The  uv  and  vw  stress  minimum  locations  are  in  close  proximity  to  one  another  and  they 

are  located  at  the  uplifting  side  of  the  outer  layer  vortical  flow  (Fig.  10).  The  Pw  stress  changes 
sign  as  the  wall  is  approached  and  shows  a  local  maximum  near  the  wall  and  a  separate  local 

maximum  closer  to  the  wing.  The  vw  stress  magnitude  is  slightly  increased  compared  to  Station 
5  values. 


CHAPTER  III.  PROBABILITY  DENSITY  FUNCTIONS  OF  u,  v, 
w  FLUCTUATIONS  AT  THE  STATIONS  SURROUNDING 

THE  VORTEX  CORE 


Inverse-instantaneous-velocity-magnitude  bias-corrected  probability  density  functions  (pdf) 
of  the  u,  V,  w  fiuctuating  velocities  are  presented  in  Figures  11  to  14.  The  pdfs  were  calculated 
at  stations  s/t=-0.247,  -0.291 ,  -0.336,  and  at  s/t=-0.368.  The  selected  stations  encompass  the 
vortex  core  located  at  s/t=-0.32,  log,o(y/f)=-1.1.  The  pdfs  show  the  probability  of  finding  a 
fluctuating  veiocity  in  a  certain  velocity  range.  The  mean  value  of  a  pdf  is  zero.  The  pdfs  were 
calculated  at  selected  stations  and  at  all  y  locations  and  equal  probability  contours  were 
generated.  The  probability  values  above  0.01  are  not  included  in  the  contour  plots  to  increase  the 
resolution  of  the  plots.  Probability  at  large  positive  or  negative  fluctuating  velocity  ranges  at  fixed 
y  location  indicates  that  the  pdfs  are  skewed  towards  positive  or  negative  side.  The  most  probable 
velocity  range  is  the  range  where  most  of  the  fluctuations  are  clustered.  The  width  of  the  pdf  of 
a  certain  fluctuating  velocity  at  a  fixed  y  location  is  an  indicator  of  the  magnitude  of  the  normal 
stress  computed  with  that  fluctuating  velocity.  The  wider  the  pdf  the  larger  the  normal  stress. 

IIM.  s/t=-0.247 

The  u  fluctuating  velocity  pdfs  show  that  the  large  negative  fluctuations  occur  in  the 
logio(y/0=-2.7  to  -2.3  range.  Most  probable  velocity  range  is  larger  than  zero  (Fig.  11a).  Very 
near  the  wall  the  most  probable  range  is  less  than  zero.  Away  from  the  wall  around  logi5(y/f)=- 
1.1  the  large  positive  and  negative  fluctuations  are  observed. 

The  V  fluctuating  velocities  at  this  station  show  both  large  positive  and  negative 
fluctuations  in  the  log,o(y/f)=-1.5  to  -1.0  range  with  larger  positive  values  (Fig.  11b). 

The  w  fluctuations  in  the  log,o(y/f)«-2.7  to  -2.0  range  show  equal  probability  velocity 
ranges  separated  with  lower  probability  velocity  ranges  (Fig.  11c).  The  positive  seldom  large  w 
fluctuations  skew  the  pdfs  towards  positive  side  and  the  most  probable  velocity  range  is  negative. 
This  indicates  that  instantaneous  velocity  prefers  two  states  or  it  indicates  the  meandering  of  the 
flow  between  two  states.  The  pdfs  around  logio(y/f)=-1-1  also  show  that  large  positive  w 
fluctuations  occur  in  this  region  and  the  most  probable  velocity  ranges  are  negative. 

III-2.  s/ts-0.291 


The  u  fluctuation  pdfs  in  the  log,o(y/f)=-2.7  to  -2.4  range  show  that  large  negative  u 


fluctuations  skew  the  pdfs  to  the  negative  side  (Fig.  12a).  Around  logi(,(y/f)=-1.1  both  large 
positive  and  negative  fluctuations  are  observed. 

The  V  fluctuations  are  spread  across  a  large  velocity  range  in  the  log,j,(y/f)=-1.5  to  -0.8 
range  (Fig.  12b).  The  positive  v  fluctuations  are  larger  than  the  negative  v  fluctuation  magnitudes. 
The  most  probable  values  in  this  range  are  negative.  Large  positive  v  fluctuations  indicate  seldom 
upwards  motion  of  the  instantaneous  flow  with  respect  to  the  mean  flow  which  might  be  related 
to  the  non-periodic  stretching  and  re-accumulation  of  the  vorticity  at  the  wall/wing  junction  at  the 
nose  of  the  wing. 

The  w  pdfs  at  this  station  show  clear  bimodal  pdfs  near  the  wall  in  the  log,o(y/ 1)~-2.7  to  - 
1.6  range  (Fig.  12c).  The  large  positive  and  negative  fluctuations  in  this  range  are  not 
accompanied  by  large  v  fluctuations.  This  indicates  that  the  meandering  of  the  vortex  between 
two  preferred  states  demonstrates  itself  with  large  variations  of  the  u  and  w  fluctuations.  The  w 
fluctuations  around  log,(,(y/f)=-1.1  are  in  a  wider  velocity  range  than  that  are  observed  at  s/t=- 
0.247.  The  pdfs  are  not  bimodal.  This  indicates  that  the  instantaneous  flow  below  the  vortex  core 
appears  to  have  preferred  directions.  The  outer  layer  vortex  appears  to  be  not  affected  by  the 
large  bimodality  of  the  near  wall  flow. 

Iil>3.  s/t=-0.336 

The  u  fluctuating  velocity  pdfs  at  this  station  also  show  large  positive  and  negative 
fluctuations  in  the  logi(j(y/f)=-2.8  to  -2.0  range  and  also  in  the  logi(,(y/f)=-1.4  to  -0.8  range 
similar  to  the  other  station  pdfs  discussed  (Fig.  13a).  Weak  bimodality  of  the  pdf  at  log^o(y/f)=*2.4 
is  an  indication  of  the  meandering  of  the  flow  in  the  y-s  plane  around  this  region. 

The  V  fluctuation  pdf  positive  values  at  logio(y/f)«-1.5  to  -0.8  range  are  smaller  in 
magnitude  than  the  values  obtained  at  station  s/t=-0.291  in  the  same  range  (Fig.  1 3b).  Positive 
V  fluctuation  values  are  larger  than  the  negative  v  fluctuation  magnitudes.  These  two  stations  are 
both  in  the  vicinity  of  the  vortex  core  and  V/U,e,  values  at  log,o(y/f)®-1,1  for  both  stations  are 
close  to  zero.  The  non-periodic  stretching  of  the  vortex  to  a  narrower  cork-screw  motion  along  the 
wall  and  wing  junction  contour  is  believed  to  cause  the  large  negative  v  fluctuations.  The  following 
expansion  of  the  vortex  into  a  larger  domain  after  the  re-accumulation  of  the  vorticity  around  the 
nose  of  the  wing  is  believed  to  cause  the  large  positive  v  fluctuations. 

The  w  pdf  bimodality  at  this  station  is  confined  to  logio(y/f)=-2.2  to  -1 .3  region  indicating 
the  extent  of  the  bimodal  flow  (Fig.  13c).  The  bimodal  region  moves  away  from  the  wall. 

ill-4.  s/t=-0.368 

The  bimodal  region  further  moves  away  from  the  wall  and  is  confined  to  log,o(y/f)=-1.7 
to  -1.3  region  (Figures  14a,  14b,  14c).  The  u  and  v  pdfs  are  similar  to  the  other  station 
distributions  with  reduced  magnitudes  of  the  fluctuations. 


CHAPTER  IV.  TRIPLE  PRODUCTS 


The  signs  of  the  Reynolds  averaged  triple  products  indicate  the  direction  of  the  average 
transfer  of  the  u^,  v^,  momentum  along  the  coordinate  axes.  The  sign  of  a  triple  product 

depends  on  the  instantaneous  flow  characteristics.  As  an  example  the  conditions  under  which  IheTFv 
is  positive  are  listed  below: 

a)  The  IFv  is  positive  when  most  of  the  large  u  fluctuations  occur  together  with  average 
magnitude  positive  v  fluctuations. 

b)  The  TFv  is  positive  when  most  of  the  large  positive  v  fluctuations  occur  together  with 
average  magnitude  u  fluctuations. 

c)  The  IFv  is  positive  when  large  u  and  large  positive  v  fluctuations  seldom  occur  together 
and  outweigh  the  seldom  occurring  large  u  and  large  negative  v  fluctuations. 

d)  Small  u  and  small  positive  v  fluctuations  outweigh  the  small  u  and  small  negative  v 
fluctuations. 

The  following  discussion  on  the  triple  products  is  based  on  the  similar  line  of  thought 
applied  to  the  other  triple  products. 

The  tPv  in  close  vicinity  of  the  vortex  core  is  large  and  positive  with  a  local  maximum 
located  at  the  uplifting  side  of  the  large  vortical  flow  (Fig.  15).  The  minimum  is  below  the  vortex 

core  at  logio(y/f)*-1.4.  Both  at  these  locations  the  values  are  also  large.  Within  the  vicinity 

of  the  vortex  core  the  tF  value  is  on  the  order  of  half  of  the  v^.  Another  local  maximum  is 

located  near  the  wall  below  the  minimum  IFv  region  at  log,o(y/f)=-2.1  and  extends  towards  the 

wing.  In  this  region  values  are  small  but  IF  values  are  large.  The  large  positive  v  fluctuations 
occur  with  average  magnitude  u  fluctuations  in  the  vicinity  of  the  vortex  core  except  right  below 
the  vortex  core.  Right  below  the  vortex  core  large  negative  v  fluctuations  occur  with  the  average 
magnitude  u  fluctuations.  Nearer  the  wall  large  u  fluctuations  occur  together  with  average 
magnitude  positive  v  fluctuations.  The  flow  that  is  uprising  encounters  the  lateral  low  turbulence 
free-stream  flow  and  results  in  positive  v  fluctuations.  On  the  other  hand  the  mean  velocity  on  the 
downward  side  directed  towards  the  wall  is  low  turbulent  free-stream  flow.  The  positive  v 
fluctuations  in  that  region  may  be  due  to  the  free-stream  velocity  that  is  aperiodically  flowing  in 
the  plane  parallel  to  the  wall  thus  causing  a  positive  v  fluctuation  on  the  average  mean  V  velocity. 
Another  explanation  may  be  the  meandering  of  the  vortex  core  in  the  plane  of  the  measurements, 
but  this  motion  is  not  very  clear  from  the  histograms  of  the  w  fluctuations.  The  minimum  right 
below  the  vortex  core  is  believed  to  be  due  to  the  fluctuations  caused  by  the  instantaneous  vortex 
structure  that  is  directed  in  to  this  region.  The  maximum  near  the  wall  is  below  y*»60  ( log,(,(y/  f)=- 


2.1).  This  region  is  where  is  low  and  it  is  close  to  the  IP  maximum  location. 

The  maximum  TPw  is  located  close  to  the  regions  where  IP  values  show  a  peak  (Fig.  1 6). 
The  transfer  of  u^  momentum  in  the  Reynolds  average  sense  is  with  positive  w  velocity.  Negative 
peak  values  are  observed  at  the  wing/wall  junction  and  at  the  uplifting  side  of  the  vortical  flow. 
The  u  fluctuations  and  positive  w  fluctuations  correlate  well  in  this  region.  The  structures  persist 
in  their  original  direction  and  directionally  lag  the  mean  flow  as  tabulated  in  Table  1.  in  Table 
1  shows  the  x  wall  coordinate  direction  with  respect  to  the  22.4°  line. 

Large  positive  TP  values  are  located  very  near  the  wall  and  the  large  negative  TP  values 
are  located  above  these  regions  (Fig.  17).  The  minimum  value  magnitude  is  twice  as  large  as  the 
maximum  value  magnitude.  The  positive  value  Indicates  an  occasional  higher  speed  fluid  reaching 
very  close  to  the  wall  resulting  in  sweeps  (+u,  -v).  The  low  momentum  fluid  lifted  up  by  these 

sweeps  result  in  slower  speed  fluctuations  causing  ejections  (-u,  +v).  The  positive  and  negative  IP 

magnitudes  are  close  to  one  another  at  the  wing  and  wall  junction.  High  TP  gradient  near  the 
wing-wall  junction  indicates  high  diffusion  rates  in  this  region. 

The  TPv,  TP  and  TPw  triple  products  together  explain  the  u^  momentum  transfer  in  the 

average  (Figures  15,  16  and  17).  The  TP  triple  product  indicates  that  the  near  wall  TP  transfer 
is  with  a  +u  transport  velocity  and  away  from  the  wall  it  is  with  a  -u  longitudinal  transport  velocity. 

The  TPv  triple  product  indicates  that  the  momentum  transfer  of  TP  except  near  the  wall  and 
below  the  vortex  core  is  with  +v  transport  velocity.  In  the  Reynolds  average  sense  the 
momentum  is  carried  towards  the  wall  in  the  near  wall  region  and  in  the  region  below  the  vortex 
core  but  carried  away  from  the  wall  within  the  vicinity  of  the  vortex  core  and  the  locations  where 

ejection  motions  are  dominant.  The  TPw  triple  product  shows  that  the  TP  momentum  is  carried 
towards  the  wing  at  the  locations  where  ejections  are  dominant.  These  iocations  are  in  the  vicinity 

of  the  vortex  core  and  below  the  vortex  core.  The  TP  momentum  is  carried  away  from  the  wing 
at  locations  very  near  the  wing-wall  intersection  and  at  the  uplifting  side  of  the  large  outer  vortical 
structure. 

The  negative  u  fluctuations  and  v  fluctuations  appear  to  be  well  correlated  in  the  region 

surrounding  the  vortex  core  as  indicated  by  the  Tjp  (Fig.  18).  The  'P  is  also  large  in  this  region. 
The  region  extends  close  to  the  wall. 

The  triple  product  TP  shows  that  in  the  vicinity  of  the  vortex  core  positive  v  fluctuations 
are  dominant  (Fig.  19).  The  maximum  value  is  obtained  at  the  downwash  side  of  the  vortical 

motion.  Large  "P  may  be  due  to  intermittent  nature  of  the  vortex  and  meandering  of  the  vortex 
location.  The  instantaneous  flow  parallel  to  the  wall  results  in  a  positive  v  fluctuation  on  the  mean 
velocity  when  the  vortex  is  instantaneously  stretched  and  confined  to  a  smaller  area. 

The  negative  w  fluctuations  and  v  fluctuations  are  well  correlated  in  the  vicinity  of  the 

vortex  core  as  shown  in  the  Pw  figure  (Fig.  20),  especially  below  the  vortex  core  and  on  the 
uplifting  side  of  the  vortex.  At  the  downward  motion  section  of  the  vortex  the  +w  fluctuations  are 


well  correlated  with  the  v  fluctuations.  The  minimum  value  magnitude  of  the  triple  product  is  twice 
as  large  as  the  maximum  value.  The  flow  below  the  vortex  core  away  from  the  wing  is  opposed 
by  an  adverse  pressure  gradient  at  s/t=-0.3. 

The  momentum  transfer  in  X22.4  direction  on  the  uplifting  side  of  the  large  vortical  flow 
is  with  a  negative  u  transport  velocity.  In  the  same  region.  Lateral  and  vertical  momentum  transfer 

is  both  away  from  the  wing  and  the  wall.  However  on  the  downward  motion  side  of  the  vortex  the'F 
momentum  is  transferred  towards  the  wing  and  longitudinal  transfer  is  with  the  negative  u 

transport  velocity.  The  vertical  momentum  transfer  is  with  a  larger  transport  velocity  than  in 
the  uplifting  region. 

The  triple  product  shows  that  in  the  vicinity  of  the  vortex  core  and  at  the  uplifting  side 
of  the  vortex  core  negative  u  fluctuations  are  more  correlated  with  the  w  fluctuations  (Fig.  21). 
Below  the  vortex  core  closer  to  the  wall  high  positive  values  are  surrounded  by  large  negative 

values.  Near  the  wall  Uw^  values  at  locations  As/t=-0.4  and  closer  to  the  wing  near  wall  uw^ 

values  are  positive.  In  this  region  sweep  motion  (+u,  -v)  is  dominant.  Near  the  separation  line  the  uw^ 
triple  product  is  small  and  positive. 

The  w  fluctuations  are  more  correlated  with  the  positive  v  fluctuations.  The  vw^  triple 
product  negative  values  are  an  order  of  magnitude  smaller  than  the  positive  maximum  value  (Fig. 
22).  High  positive  values  are  located  around  the  uplifting  region  of  the  vortex.  The  small  negative 
values  are  observed  in  the  near  wall  region  below  the  vortex  core  and  in  the  outer  edge  of  the 
downward  side  of  the  vortex. 

The  vP  triple  product  maximum  and  the  minimum  peak  locations  are  located  below  the 

vortex  core  (Fig.  23).  Large  negative  W  are  located  at  the  uplifting  side  of  the  vortex  and  large 

positive  values  are  observed  at  the  downward  side  of  the  vortex.  Positive  values  show  the 
seldom  but  strong  positive  w  fluctuations  are  dominant  and  the  w^  momentum  is  carried  with  a 
positive  lateral  transport  velocity  in  the  positive  s  direction.  The  elevated  values  are  believed  to 
be  due  to  the  meandering  of  the  vortex  core. 

The  triple  products  UvP,  vw^  and  wP  together  describe  the  transfer  of  the  "vP  momentum 

(Figures  21 , 22  and  23).  The  Uw^  shows  that  while  the  W  momentum  is  transferred  with  positive 

u  transport  velocity  in  a  column  below  the  vortex  core  TF  momentum  is  transported  with  negative 

u  velocity  around  this  column  and  in  the  vicinity  of  the  vortex  core.  The  'vvP  values  show  that  in 
the  vicinity  of  the  vortex  core  and  at  the  uplifting  side  of  the  large  outer  vertical  motion  momentum 

is  carried  away  from  the  wall.  This  region  is  also  where  negative  u  fluctuations  are  dominant.  The'vF 

momentum  carried  away  from  the  wall  is  transported  with  negative  u  transport  velocity.  ThelF 
triple  product  on  the  other  hand  shows  that  the  momentum  is  transferred  both  towards  the  wing 
and  away  from  the  wing  due  to  the  meandering  of  the  instantaneous  vortex  structure. 

The  triple  product  uvw  (Fig.  24)  is  positive  when  seldom  but  positive  fluctuations  occur 
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simultaneously  such  that  the  average  value  is  tipped  towards  the  positive  or  when  two  of  the 
fluctuating  velocities  have  seldom  but  strong  negative  fluctuations  provided  that  the  third 

fluctuating  velocity  is  simultaneously  positive  and  strong.  The  uvw  is  positive  in  the  uplifting  side 

of  the  vortex  and  largest  negative  values  occur  beneath  the  vortex.  Large  positive  Uvw  values 
are  also  located  near  the  wall.  Large  negative  values  indicate  that  large  seldom  negative 

fluctuations  or  two  positive  and  one  negative  strong  fluctuations  occur  together.  The  TF  triple 
product  is  positive  at  the  regions  where  uvw  is  large. 


CHAPTER  V.  DERIVED  QUANTITIES 


V-1.  FLOW  ANGLE 

The  flow  angle  is  defined  with  respect  to  the  X22.4  axis.  A  positive  flow  angle  shows  a  fluid 
motion  towards  the  wing  (Fig.  25).  The  maximum  flow  angle  is  at  the  boundary  layer  edge  of 
station  5.  The  large  vortex  structure  results  in  positive  flow  angles  above  the  vortex  core  and 
negative  flow  angles  below  the  vortex  core.  Near  the  wall/wing  junction,  positive  flow  angles  are 
observed.  The  minimum  value  is  located  between  a  vertical  line  passing  through  the  vortex  core 
and  the  wing  near  the  wall.  The  flow  angles  are  of  opposite  signs  on  each  side  of  the  separation 
zone,  indicating  that  the  flow  is  merging  towards  a  region  from  both  sides. 

V-2.  FLOW  GRADIENT  ANGLE 

Flow  gradient  angle  (FGA)  values  show  large  variations  between  -114°  and  90°  with 
respect  to  the  X22.4  line  (Fig.  26).  A  minimum  value  is  obtained  away  from  the  vortex  core,  wing 
and  the  wall,  and  it  is  located  below  the  vortex  at  s/t=-0.35,  logio(y/t)=-1 .6.  The  FGA  shows  large 
gradients  between  logiQ(y/t)*-1 .6  and  -2.2.  Positive  FGA  values  are  located  in  the  zone  between 
the  vortex  and  the  wing.  Sharp  variations  are  also  observed  near  the  wing. 

V-3.  SHEAR  STRESS  ANGLE 

The  shear  stress  angle  (SSA)  has  large  positive  values  beneath  the  vortex  and  large 
negative  values  near  the  wall  in  the  X22.4-S  plane  of  the  measurements  (Fig.  27).  Negative  and 

positive  values  greater  than  90°  indicate  a  region  where  -  ov  shear  stress  is  negative.  A  positive 

SSA  shows  a  -  vw>0  and  negative  SSA  indicate  a  -  vw<0.  The  plot  shows  that  very  near  the  wall 

at  the  separation  location,  the  -uv  shear  stress  is  negative.  In  such  a  flow  the  SSA  lags  highly 
behind  the  FGA.  Large  variations  by  the  FGA  are  not  followed  by  the  SSA. 

V-4.  (dW/dy)  t/U^ 

The  vorticity  component  in  the  X22.4  direction  is  (dW/dy-BV/ds)  (Fig.  28).  In  the  figure  the 
larger  of  the  components  dW/dy  is  presented.  The  figure  shows  that  large  outer  layer  vortex  has 
positive  vorticity  with  a  small  magnitude.  Vorticity  is  the  rotation  of  the  fluid  particles  along  an  axis 
and  does  not  necessarily  imply  a  large  structural  rotation.  Maximum  negative  vorticity  is  located 
very  near  the  wall  close  to  the  wing  wall  junction.  It  appears  that  the  large  outer  layer  vortex 


generates  a  near  wall  vortex  with  the  opposite  sense  of  rotation.  At  the  wall/wing  junction  the 
vorticity  values  reach  to  their  positive  peak  with  the  same  sense  of  vorticity  of  the  outer  region 
vortex  and  are  caused  by  the  large  negative  vorticity  region  near  the  wall.  The  vorticity  values 
around  the  rotational  separation  zone  are  negative. 

V-5.  TKE/U„," 

The  turbulent  kinetic  energy  is  defined  as: 

2 

Figure  29  shows  that  the  TKE  is  very  large  in  the  region  below  the  core.  In  this  complex 
flow  the  TKE  levels  are  very  large  away  from  the  wall.  In  a  3-D  turbulent  boundary  layer  flow  TKE 

maximum  value  is  usually  at  y*«25.  Large  TKE  values  are  due  to  large  W  values.  In  the  vicinity 
of  the  rotational  separation  region  TKE  values  are  higher  than  at  Station  5. 

V-6. 

V^2  is  defined  as  the  transport  velocity  of  the  turbulent  kinetic  energy: 

The  ratio  defines  a  Reynolds  average  sense  perpendicular  velocity  component  that  the 
TKE  is  transported  (Fig.  30).  The  velocity  is  observed  to  be  positive  in  general.  Large  positive 
values  are  located  in  the  vicinity  of  the  vortex  core  both  on  the  uplifting  and  downwash  side  of 
the  vortex.  The  TKE  is  transported  away  from  the  wall  in  these  regions.  Large  negative  values 
are  located  at  the  downward  side  of  the  vortex.  This  indicates  that  the  TKE  generation  is 
maximum  away  but  close  to  the  wall  and  the  TKE  is  transported  in  both  directions  around  this 
region. 

V-7.  uvfiF 

The  structural  parameter  shows  that  the  contours  define  plateaus  near  the  wall  region 
below  log,o(y/ 1)~  -2  indicating  that  the  IP  is  the  correct  scaling  parameter  for  the  uv  near  the  wall 
even  in  highly  complicated  flows  (Fig.  31).  This  information  can  be  used  to  deduce  theluv 

information  near  the  wall  from  measured  TF  data.  However  the  coefficient  values  are  not  high 
indicating  that  the  u  and  v  fluctuating  velocities  are  not  very  well  correlated. 

V-8.  SKEWNESS  OF  u 

Skewness  of  the  u  fluctuating  velocity  is  defined  as: 
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(7F)^ 

The  quantity  shows  the  skewing  of  the  probability  density  functions  due  to  seldom  but 
strong  fluctuations  (Fig.  32).  The  skewness  of  a  Gaussian  pdf  is  zero.  The  is  larger  than  zero 
near  the  wall  indicating  that  the  seldom  but  strong  positive  u  fluctuations  occur  in  this  region. 
These  are  due  to  the  sweeps  which  are  directed  towards  the  wall  originating  from  a  region  of  high 
momentum  away  from  the  wall.  Large  negative  values  near  the  wing  are  due  to  the  large  vortical 
flow  generated  near  the  wing  reorienting  itself  due  to  the  presence  of  the  wing  and  the  wall.  In 
the  vicinity  of  the  vortex  core  negative  u  fluctuations  dominate  the  flow. 

V>9.  SKEWNESS  OF  v 

Large  positive  skewness  values  are  located  above  the  vortex  core  with  the  maximum 
values  located  at  the  uplifting  side  of  the  vortex  (Fig.  33).  Large  negative  values  are  located  at 
the  downward  side  of  the  vortex.  This  shows  that  the  seldom  but  strong  v  fluctuations  occur  at 
the  uplifting  side  and  above  the  vortex  core.  The  flow  on  the  downwash  side  is  occasionaliy  with 
a  fast  instantaneous  velocity.  In  general  the  is  larger  than  zero.  Negative  values  are  located 
near  the  wall  and  away  from  the  vortex  core  and  in  the  vicinity  of  the  rotational  separation  region. 
Negative  S^  in  this  region  shows  that  the  flow  occasionally  moves  upwards  with  a  fast 
instantaneous  velocity. 

V-10.  SKEWNESS  OF  w 


Skewness  of  w  fluctuations  S^  shows  that  large  positive  values  are  located  below  the 
vortex  core  all  the  way  to  the  wall  (Fig.  34).  On  the  contrary  to  the  intuitive  reasoning  the 
skewness  in  this  region  is  positive.  The  positive  values  indicate  that  occasional  strong  +w 
fluctuations  occur  below  the  vortex  core  region.  This  may  be  due  to  chaotic  meandering  of  the 
vortex  zone.  Between  this  zone  and  the  wing  negative  S„  indicates  occasional  strong  -w 
fluctuations.  Large  positive  values  show  that  the  occasional  large  positive  and  small  negative 
fluctuations  are  caused  by  the  flow  entrained  by  the  vortex.  The  negative  values  located  under 
the  vortex  core  near  the  wall  is  also  believed  to  be  due  to  large  chaotic  meandering  of  the  vortex 
structure.  Negative  values  surround  the  separation  region  and  near  the  wall  to  indicate  that 
the  redirected  flow  generates  seldom  but  strong  -w  fluctuations.  Near  wall  values  on  the  left-hand 
side  of  the  separation  are  small  but  positive. 

V-11.  SHEAR  STRESS  MAGNITUDE 

Shear  stress  magnitude  is  highly  elevated  in  the  return  and  uplifting  side  of  the  vortex  (Fig. 
35).  The  fluctuations  become  more  correlated  in  these  regions  since  -w  and  +v  fluctuations 
describe  the  instantaneous  motion  of  the  fluid  redirected  by  the  vortex  in  this  region.  The  -u  and 
+v  fluctuations  are  also  large  in  these  regions.  The  separation  region  is  observed  to  have  the 


lowest  shear  stress  magnitude  but  larger  than  Station  5  values. 


V-12.  TOWNSEND’S  STRUCTURAL  PARAMETER,  A, 

Townsend’s  structural  parameter  is  defined  as: 

IF  +'iF + mF 

The  A,  structural  parameter  increases  from  very  small  values  near  the  wall  to  large  values 
at  the  outer  region  (Fig.  36).  The  downward  side  of  the  vortex  core  region  is  observed  to  have 
small  values  due  to  small  shear  stresses  in  this  region  since  the  free-stream  fluctuations  are  close 
to  the  isotropic  distribution.  However  large  values  are  observed  on  the  uplifting  side  of  the  vortex. 

V-13.  STRUCTURAL  PARAMETER 

■F 

The  values  are  close  to  1  in  most  part  of  the  layers.  This  value  is  the  half  of  the  isotropic 

distribution  value  (Fig.  37).  This  shows  that  F  fluctuations  are  in  general  much  higher  than  the 
isotropic  value  away  from  the  wall  which  is  due  to  large  vortical  motion  energizing  the  whole 
domain  studied.  Near  the  wall  region  especially  close  to  the  wall/wing  junction  the  values  are 
much  larger  than  2  which  is  mainly  due  to  very  high  u  fluctuations  dominating  the  region. 

V-14.  ROTTA’S  T 

Rotta’s  T  anisotropy  constant  is  defined  as: 

T^tan(SSA-FA) 

tan(FGA-FA) 

The  eddy  viscosity  ratio  T  is  relative  to  the  local  flow  angle  and  is  co-ordinate  system 
independent  about  the  y-axis.  T  is  not  a  constant  at  any  y  for  these  data  according  to  Fig.  17.  T 
is  about  constant  in  the  outer  layer  with  T«0.4  to  0.5.  Wherever  FGA=FA,  the  data  show  T 
approaching  large  negative  or  positive  values.  Below  log, ^(y/f)  of  about  -1 .5,  this  is  a  poor  eddy- 
viscosity  ratio  model. 


CHAPTER  VI.  VARIATION  IN  FLOW  PHYSICS  AT 
SELECTED  LOCATIONS 


VI-1.  SEPARATION  LOCATION 


In  this  section  the  profiles  of  flow  quantities  measured  at  stations  5,  A,  B  and  C  are  used 
to  describe  the  flow  around  a  three-dimensional  separation  region.  Mean  secondary  flow 
streamlines  in  the  plane  of  the  measurements  around  s/t=-0.58  show  a  region  where  the 
secondary  streamlines  merge.  The  streamlines  diverging  away  from  the  wall  join  to  this  rotational 
region  which  indicates  separation  of  the  flow  from  the  floor  around  s/t=-0.54.  The  vorticity  of  the 
region  is  opposite  to  that  of  the  large  outer  layer  vortex. 

Within  the  surrounding  of  the  rotational  separation  region  the  mean  velocity  magnitude  is 
small  since  the  region  is  close  to  the  wall.  Compared  to  station  5  skin  friction  velocity  u,=1.15 
m/sec  on  the  wall,  u^  decreases  to  1 .094  m/sec  at  station  A,  to  1 .002  m/sec  at  station  B,  and  to 
1 .048  m/sec  at  station  C.  The  W  mean  velocity  minima  values  decrease  at  stations  A,  B  and  C, 
and  the  minima  locations  shift  up  in  the  layers.  Station  A  appears  to  be  the  station  closest  to 
where  the  rotational  separation  region  away  from  the  wall  occurs.  The  oil  flow  visualization  on  the 
floor  shows  that  station  C  is  closer  to  the  separation  on  the  wall  where  the  streamlines  start  to 
converge. 

The  TF  normal  stress  (Fig.  5)  at  logio(y/f)=-2.5  increases  14%  by  station  C  with  a  larger 


increase  nearer  the  wall.  The  wF  normal  stress  (Fig.  7)  steeply  increases  near  the  wall 
proceeding  from  station  5  to  station  C.  The  monotonic  increase  show  a  100%  increase  at  station 

C  compared  to  that  at  station  5  at  log,o(y/0=-1-5.  The  "F  peak  value  is  also  increased 


proceeding  towards  the  wing  (Fig.  6).  The  peak  value  in  the  outer  layer  is  doubled.  The  uv  (Fig. 

8)  and  vw  (Fig.  10)  stress  minimum  and  maximum  value  magnitudes  increase  within  this  zone. 
The  u  and  w  and  the  v  and  w  fluctuations  are  more  correlated  within  the  separation  region 
compared  to  Station  5. 

The  transport  of  the  IP  momentum  both  with  positive  and  negative  u,  v,  w  fluctuations  is 


increased.  The  maxima  and  minima  value  magnitudes  of  TPv  (Fig.  15)  ,  TFw  (Fig.  16)  andlF 
(Fig.  17)  increase.  The  increase  in  magnitudes  at  the  same  y/t  locations  away  from  wall  depicts 


larger  diffusion  for  the  transport  of  the  IP 


normal  stress  (TD=- 


dXf 


(lPUi),  largest  term 


TD=- 


The  "P  momentum  is  transported  increasingly  with  negative  u  and  w  transport 


velocities  and  transported  increasingly  with  positive  v  transport  velocity.  Diffusion  of  is 

increased  (TD=  largest  term=  -^(  v^)  (Fig  19)).  The  uw^  minima  and  maxima  vaiues 

OX/  dy 

increase,  ~\fvF  (Fig.  22)  is  positive  and  maxima  vaiues  increase,  become  increasingly 

negative.  These  observations  show  that  TP  is  increasingly  transported  away  from  the  wall  and 

the  wing.  The  longitudinal  transport  of  TP  is  increased  both  with  positive  and  negative  u 
fluctuations. 

VI-2.  STATIONS  SURROUNDING  THE  VORTEX  CORE 

The  stations  discussed  in  this  section  are  TB,  TC,  R  and  S  which  encompass  the  vortex 
core  of  the  outer  layer  vortex.  Even  though  the  streamlines  indicate  a  rotating  structure  in  the 
outer  layer,  the  vorticity  of  the  flow  in  this  region  is  actually  small.  However,  the  circulation  is  large 
since  it  covers  a  large  area.  The  near  wall  flow  affected  by  this  large  structure  shows  a  large 
vorticity  region  in  a  very  narrow  region  near  the  wall  with  the  opposite  sense  of  rotation  to  that 
of  the  large  vortex.  The  vortex  core  is  located  at  log,o(y/f)=-1 .1 .  Below  the  vortex  core,  U  mean 

velocity  profiles  (Fig.  2)  show  regions  where  9iy/3y  is  zero  or  even  negative  which  indicates  the 
largest  gradient  term  using  boundary  iayer  approximations  is  nearly  zero  for  such  a  flow.  The 
maximum  V  velocities  (Fig.  3)  are  observed  near  the  vortex  core  and  W  mean  velocities  (Fig.  4) 
change  sign  above  the  vortex  core.  The  regions  where  dU/dy=0  (between  logio(y/f)=-2  to  -1 .25) 

shows  that  the  production  of  tP  and  Uv  stresses  are  small.  In  these  regions  IF  normal  stress 
values  decrease  proceeding  from  station  TB  to  station  S.  The  dominant  normal  stress  near  the 

vortex  core  is  "P  with  values  larger  than  near  wall  IF  values.  The  maximum  value  of  the  TP  near 
the  wall  is  much  larger  than  the  other  normal  stresses  due  to  aperiodic  motion  of  the  near  wall 
flow  to  generate  bimodal  w  fluctuating  velocity  histograms.  The  w‘  fluctuations  reach  0.23  U,,,  at 

logio(y/0=-2  and  v‘  fluctuations  are  as  high  as  0.19  U„,  at  logio(y/f)=-1.15.  The  TP  andP 
stresses  decrease  on  both  sides  of  the  vortex  core. 

The  Uv  correlation  obtains  the  largest  values  within  the  measurement  plane  at  station  TB 
with  values  on  the  order  of  3.5  larger  than  at  station  5.  However  the  values  rapidly  decrease  by 

station  S.  The  locations  where  dUldy=0  also  corresponds  to  the  locations  where  -P.^  which 

ay 

is  the  largest  Vv  production  term  is  approximately  zero.  In  this  region  -Uv—  production  term 

ay 

becomes  important.  The  uiv  stress  (Fig.  9)  on  the  uplifting  and  on  the  downwash  side  of  the 
vortex  show  different  features.  On  the  uplifting  side  the  positive  maximum  is  located  in  between 
two  minima. 

Positive  Uw  correlation  denotes  the  -(pmi^  stress  slowing  down  the  flow  to  rotate  along 
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the  -y  axis  and  a  negative  uw  correlation  denotes  a  stress  slowing  down  the  flow  to  rotate  along 
the  +y  axis.  On  the  downwash  side  of  the  vortex,  the  near  wall  flow  encounters  stresses  slowing 

down  the  flow  towards  the  vortex.  The  vw  correlation  reaches  a  minimum  value  at  station  TB 
within  the  measurement  plane  and  increases  on  both  sides  of  this  station.  The  maximum  value 

is  reached  at  station  R.  Negative  vw  results  in  a  rolling  motion  along  the  +X22.4  axis;  this  stress 
is  the  magnitude  of  the  resistance  of  the  fluid  to  a  motion  along  the  +s  axis. 

The  TF  (Fig.  17)  is  positive  very  near  the  wall  at  all  stations  of  this  zone  and  above 

log,o(y/f)=-2.5  the  values  are  negative.  Within  the  same  near  wall  region  TFv  (Fig.  15)  values 
are  negative  indicating  that  sweep  type  motions  (u>0,  v<0)  are  dominant  in  this  region.  Above  the 

near  wall  region  IP  is  negative  throughout  the  layers.  Except  at  stations  on  each  side  of  the 

vortex  core  (TC  and  R)  the  TFv  is  positive  above  this  near  wall  region.  The  TF  momentum  in  this 

region  is  carried  away  from  the  wall  with  negative  u  transport  velocity.  The  TFw  (Fig.  16)  value 
becomes  positive  proceeding  from  station  TB  towards  station  S  and  the  maximum  value 

increases,  indicating  IF  momentum  transport  with  a  positive  w  transport  velocity  nearer  the  wing 
side  of  the  vortex  core. 

The  transport  of  IP  is  with  a  positive  w  transport  velocity  as  depicted  by  v^w  (Fig.  20). 
Near  the  vortex  core  at  stations  TB,  TC  and  R  it  is  with  negative  w  transport  velocity.  Transport 

velocity  is  practically  zero  above  the  vortex  core.  The  transport  of  IP  similar  to  that  of  the  IP  is 

with  negative  u  transport  velocity  as  shown  by  UP.  At  the  downwash  side  the  value  becomes 

positive  below  the  vortex  core.  The  IP  indicates  that  the  'P  transport  is  away  from  the  wall.  The  17^ 
momentum  is  transported  away  from  the  wall  and  wing  with  a  negative  u  transport  velocity. 

The  (Fig.  23)  shows  that  transport  of  vF  is  with  negative  w  transport  velocity  on  the 
uplifting  side  of  the  vortex  and  with  positive  w  transport  velocity  on  the  downwash  side  of  the 

vortex.  Transport  velocities  at  different  stations  change  sign  many  times.  The  IP  momentum  is 
also  transported  away  from  the  wall  in  a  large  part  of  the  layer  similar  to  that  of  the  and  yF 

momentum  transport.  The  v  transport  velocity  of  IP  decreases  below  the  vortex  core.  The 

transport  of  W  in  the  longitudinal  direction  is  with  a  negative  u  transport  velocity  at  station  TB 
however  the  value  becomes  positive  and  large  positive  at  stations  TC  and  R  and  at  station  S  it 

is  with  a  large  negative  u  transport  velocity.  The  IP  momentum  transport  is  strongly  affected  by 
the  vortex  motion  and  responds  rather  quickly  in  such  a  short  spatial  variation. 

VI-3.  NEAR  THE  WING/WALL  JUNCTION 

In  this  section  the  flow  physics  near  the  wing/wall  junction  is  described.  The  measurement 
stations  used  in  this  section  are  stations  T  through  X.  The  mean  U  velocity  accelerates  as  the 
wing  is  approached.  Near  the  wall,  the  velocity  gradient  is  increased.  The  region  with  constant 


U  velocity  extends  a  longer  range  than  around  the  vortex  core.  The  V  vertical  velocity  decreases 
with  the  minimum  at  station  X.  The  W  velocity  minima  values  increase  and  the  maxima  values 
decrease  as  the  wing  is  approached.  At  station  X  closest  to  the  wing  the  W  component  changes 
sign  near  the  wall. 

The  near  wall  IF  normal  stress  values  continuously  increase  as  the  wing  is  approached. 
However  the  values  above  log,o(y/f)=-2  continuously  decrease  till  station  W  to  reach  values 

below  the  TF  normal  stress  obtained  at  station  5.  At  station  X  the  values  abruptly  increase.  The 

and  "vF  normal  stresses  are  reduced  till  station  W  throughout  the  layers.  The  values  of  the  normal 
stresses  at  station  W  are  relatively  higher  than  at  station  5.  At  station  X  values  increase.  The 
decrease  of  the  normal  stresses  within  this  zone  is  due  to  the  strong  acceleration  of  the  flow  and 
also  due  to  the  free  stream  low  turbulence  fluid  carried  close  to  the  wall  by  the  large  outer  region 
vortex.  However,  at  station  X  the  effect  of  the  wing/wall  junction  is  felt.  The  turbulent  kinetic 
energy  of  the  flow  away  from  the  wall  is  lower  at  locations  where  all  normal  stresses  are  reduced. 

The  Vv  stress  values  at  stations  T  and  U  are  comparable  to  the  values  at  station  5.  At 
station  W  the  W  stress  magnitude  is  decreased  above  log,o(y/f)=-2.25  with  values  close  to  zero. 
Near  wall  values  are  as  large  as  the  values  at  station  5.  The  vw  stress  resembles  a  wave  motion 
with  increasing  and  decreasing  values  above  and  below  zero  values.  The  vw  stress  magnitude 

decreases  as  the  wing  is  approached.  At  station  X  the  vw  stress  once  more  increases.  Theuw 
stress  also  show  a  wavy  character  with  the  maximum  magnitude  value  located  nearer  the  wall. 
The  reduction  of  this  magnitude  at  station  U  compared  to  that  of  T  is  followed  by  an  increase 

which  results  in  values  an  order  of  magnitude  larger  at  station  W  than  that  at  station  5.  The TTw 

stress  changes  sign  near  the  wall  at  station  X.  The  increase  of  the  maximum  -Uw  causes 
increase  resistance  to  a  rotation  along  the  +y  axis. 

As  indicated  by  the  TFv,  IFw  and  TF  triple  products,  the  near  wall  Reynolds  averaged 

transport  of  IF  momentum  is  towards  the  wall  and  away  from  the  wing  with  positive  longitudinal 

transport  velocity.  Away  from  the  wall  TF  transport  is  away  from  the  wall  towards  the  wing  with 
negative  longitudinal  transport  velocity.  Near  the  wall  sweep  motions  are  dominant  where  positive 
large  u  fluctuations  occur  with  large  negative  v  fluctuations.  Above  this  region  ejection  events 
(u<0,  v>0)  become  dominant.  All  three  triple  products  approach  zero  above  log,o(y/f)*-2.0  as  the 
wing  is  approached.  The  triple  products  at  station  X  do  not  follow  the  trends  of  the  profiles  of 

stations  T  through  W,  especially  above  log,(,(y/f)®-2.0.  The  TFv  is  larger  than  near  the  wall  in 
this  region,  the  TFw  becomes  negative  and  the  TF  is  negative. 

The  'F  momentum  transfer  near  the  wall  is  with  positive  longitudinal  velocity  and  above 
log,o(y/0“*2-5  it  is  with  negative  longitudinal  velocity.  The  lateral  Reynolds  average  sense 

transport  velocity  of  v^  momentum  changes  sign  within  the  layers  many  times.  At  station  X  while 


the  near  wall  momentum  is  transported  towards  the  wall  above  log,o(y/f)=-1 .5  it  is  transported 

away  from  the  wall.  The  vertical  transport  velocity  of  the  momentum  also  shows  a  wave  type 
variation  in  the  layers.  At  stations  V  and  W  the  Reynolds  averaged  transport  velocity  is  positive 
throughout  the  layers  at  station  X  near  wall  magnitude  of  this  transport  velocity  is  increased 
however  in  the  outer  region  it  becomes  negative. 

The  Tiw^,  vw^  and  TF  triple  products  describe  the  transport  of  W.  The  transport  ofTF 
near  the  wall  is  with  positive  longitudinal  velocity,  above  logio(y/f)*-2.5  it  is  with  negative 
longitudinal  velocity,  in  this  region  negative  large  occasional  u  fluctuations  tip  the  average  towards 
negative  values.  Transport  velocity  minimum  location  shifts  up  in  the  layer  at  station  X.  The 

vertical  transport  of  W  is  with  positive  transport  velocity  and  in  the  outer  region  it  is  towards  the 
wall.  The  lateral  transport  shows  wavy  character  of  the  transport  with  negative,  positive  transport 
velocities  throughout  the  layer. 


CHAPTER  VII.  CONCLUSIONS 


In  the  present  work  structural  features  of  a  complex  three-dimensional  flow  in  a  wing-wall 
junction  were  studied.  The  flow  is  complex  not  only  due  to  three-dimensionality  of  the  boundary 
layer  on  the  wall  but  also  due  to  the  large  vortical  structure  present  in  the  outer  layer  and  due  to 
the  separation  of  the  3-D  TBL  affected  by  the  adverse  lateral  pressure-gradients.  The  analysis 

is  further  challenging  due  to  the  aperiodic  motion  of  the  vortex  which  generates  large  TP  normal 
stresses. 

The  large  vortical  structure  in  the  outer  region  brings  high  speed  free-stream  fluid  very 
close  to  the  wall,  which  results  in  thinning  of  the  boundary  layer  near  the  wing.  Adverse  lateral 
pressure  gradients  result  in  a  3-D  separation  zone.  The  flow  very  near  the  wing  has  positive  flow 
angle  with  respect  to  X22.4.  The  outer  layer  vortical  flow  has  small  vorticity  along  the  +X22.4.  The 
high  vorticity  region  is  concentrated  in  a  region  near  the  wall  between  the  wing  and  the  vortex 
center  with  a  vorticity  along  the  axis.  At  the  wing/wall  junction  a  region  with  vorticity  along 
the  +X22.4  axis  is  observed. 

The  TF  normal  stress  is  higher  in  the  region  between  the  vortex  core  and  the  wall  than 

the  other  regions  of  the  flow.  The  maximum  IF  is  very  near  the  wing/wall  junction.  TheTF 

maximum  occurs  below  the  vortex  core  and  large  P  values  occur  around  the  vortex  core.  TheTP 

is  very  large  compared  to  IF  and  P  normal  stresses,  especially  in  the  region  between  the  vortex 
core  and  the  wall.  Bimodal  w  fluctuation  histograms  are  observed  at  four  measurement  stations, 
two  on  each  side  of  the  vortex  core.  The  turbulent  kinetic  energy  of  the  flow  increases  to  5  to  6 
times  to  that  of  the  station  5  value  at  the  stations  in  a  region  below  the  vortex  core. 

Both  the  uv  and  the  vw  stresses  reach  their  maximum  values  on  the  uplifting  side  of  the 

vortex  core.  The  Tiw  stress  maxima  are  located  near  the  wall  below  the  vortex  core  closer  to  the 
wing  and  at  the  wing/wall  junction.  Near  the  vortex  core  and  the  separation  region  the  flow  is 
along  the  X22.4  axis.  The  inrush  of  the  free-stream  flow  towards  the  wall/wing  junction  generates 
flow  away  from  the  wing  with  flow  angles  as  high  as  -23.5°  with  respect  to  X22.4  axis.  The  shear 
stress  angle  severely  lags  behind  the  flow  gradient  angles.  Where  there  are  severe  FGA 
variations,  the  SSA  response  in  the  outer  layer  is  not  followed  in  the  inner  regions.  Due  to  small 
dU/dy  values  observed  throughout  the  layers  the  FGA  are  very  large. 

Diffusion  processes  are  highly  altered  by  the  outer  layer  vortex  effects  on  the  inner  layer 
flow.  Large  triple  product  magnitudes  are  observed  surrounding  the  vortex  core,  below  the  vortex 
core,  within  the  region  between  the  vortex  core  and  wing  nearer  the  wall  and  near  the  wing/wall 
junction. 
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TABLE  1.  LASER-DOPPLER  VELOCIMETER  LOCATIONS  AND  FLOW  PARAMETERS. 


Stations 

X 

(mm) 

z 

(cm) 

Ue 

(m/sec) 

(m/sec) 

Pw 

(degree) 
in  22.4° 

Pfs 

(degree) 
in  22,4° 

6 

(mm) 

s/t 

(ac,/ax) 

(1/m) 
in  22.4° 

Oc,/az) 

(1/m) 
in  22.4° 

Station  5 

6.60 

-7.47 

29.54 

1.15 

2.7 

14.7 

40 

-0.6774 

-5.551 

-3.471 

A 

8.79 

-6.94 

29.74 

1.094 

2.276 

13.12 

39.6 

-0.589 

-6.402 

-5.195 

B 

9.98 

-6.64 

29.94 

1.002 

-0.647 

12.53 

39.9 

-0.545 

-7.025 

-5.994 

C 

11.2 

-6.35 

30.2 

1.048 

-1.195 

12.02 

39.3 

-0.501 

-7.701 

-6.671 

D 

12.42 

-6.06 

30.45 

1.093 

-2.409 

11.41 

39.2 

-0.456 

-8.249 

-7.983 

E 

13.61 

-5.77 

30.65 

1.132 

-3.184 

10.88 

40.3 

-0.412 

-8,782 

-9.488 

TB 

14.83 

-5.47 

31.82 

1.163 

-3.858 

9.83 

39.7 

-0.368 

-9.754 

-10.455 

TC 

15.7 

-5.26 

32.27 

1.188 

-7.074 

9.29 

40.9 

-0.336 

-10.884 

-9.813 

R 

16.92 

-4.97 

32.58 

1.297 

-15.583 

8.21 

39.2 

-0.291 

-11.654 

-10.579 

S 

18.11 

-4.67 

33.15 

1.362 

-20.741 

7.09 

40.3 

-0.247 

-12.193 

-5.293 

T 

19.33 

-4.38 

33.72 

1.413 

-18.914 

5.77 

38.9 

-0.203 

-14.54 

-4.502 

U 

20.55 

-4.09 

34.85 

1.320 

-13.539 

4.24 

39.3 

-0.159 

1  -15.129 

-5.093 

V 

21.74 

-3.79 

36.29 

1.454 

-10.943 

2.82 

38.4 

-0.114 

-13.134 

-8.909 

W 

22.96 

-3.50 

37.56 

1.656 

-10.888 

1.02 

28.6 

-0.070 

1.456 

-0.949 

X 

24.26 

-3.18 

40.42 

1.917 

10.752 

-0.17 

20.6 

-0.023 

1.472 

-0.933 

Maximum  thickness  of  the  wing  (t)  =  7.17  cm,  Nominal  reference  velocity  (U,„)  =  27.5m/sec,  Reynolds  number  based 
on  momentum  thickness  at  0.75  chord  upstream  of  the  wing  on  the  tunnel  centerline  (Re,)  =  5936.  Wing  chord  length 
=  30.48  cm. 


TABLE  2. 2 
MEAN  VEL 
TRIPLE  PF 

n:1  ODDS  ±2 
.OCITIES,  RE> 
lODUCTS. 

a  UNCERTAINTIES  OF 
^NOLDS*  STRESSES  AND 

Term 

Uncertainty 

X  10-“ 

Term 

Uncertainty 

X  10-°® 

u/u,„ 

5.64 

4.54 

v/u,„ 

1.25 

7.97 

w/u,„ 

3.7 

2.82 

0.09 

3.04 

0.0581 

9.13 

0.01017 

2.93 

WlXiJ 

0.0363 

2.77 

uwlUj 

0.0537 

7F/u„® 

14.06 

0.0276 

■p/u,„® 

2.38 

Wi\iJ 

23.47 

TABLE  3.  SOME  LENGTH  SCALES  OBTAINED  FROM  LDV  DATA  IN  22.4“  LINE  COORDINATES 


63(mm) 


4.85 


6(mm) 

5,  (mm) 

52(mm) 

40.0 

6.17 

5.26 

39.6 

6.02 

5.27 

39.9 

5.95 

5.26 

39.3 

5.81 

5.16 

39.2 

5.71 

5.09 

40.3 

5.71 

5.1 

39.7 

5.51 

4.94 

40.9 

5.37 

4.79 

39.2 

5.01 

4.42 

40.3 

4.64 

4.14 

38.9 

4.04 

3.69  : 

39.3 

3.72 

3.49  ; 

38.4 

3.46 

3.33  ; 

28.6 

2.46 

2.42  ; 

20.6 

1.82 

1.78 

5.,(mm) 

55(mm) 

56(mm) 

57(mm) 

-7.86 

0.79 

-7.07 

-1.68 

-6.9 

0.57 

-6.32 

-1.38 

-6.47 

0.49 

-5.98 

-1.29 

-6.08 

0.41 

-5.67 

-1.21 

-5.76 

0.34 

-5.42 

-1.14 

-5.35 

0.27 

-5.08 

-1.11 

-4.74 

0.17 

-4.57 

-1.02 

-4.49 

0.14 

-4.36 

-1.02 

-3.9 

0.01 

-3.89 

-1.00 

-3.22 

-0.03 

-3.26 

-0.85 

-2.53 

-0.05 

-2.59 

-0.61 

-1.78 

-0.06 

-1.84 

-0.39 

-1.20 

-0.06 

-1.25 

-0.23 

-0.31 

-0.07 

-0.38 

-0.08 

0.76 

-0.09 

0.67 

-0.07 

Ug=Velocity  magnitude  at  the  layer  edge. 

6=  y  where  71/^=0.995=  Boundary  Layer  Thickness 


6,=  r(1--ii)dy=  Streamwise  Displacement  Thickness 

*L  U 

0 

62=  (1“i - )Gfy=  Magnitude  Displacement  Thickness 

u 

0 

53=J(1  -^)^dy=  Streamwise  Momentum  Thickness 

7  IV 

64=  (“ — )c/y=  Lateral  Displacement  Thickness 

^  u 

0  e 

J"  w  u 

— (1  - — )dy=  Lateral  Momentum  Thickness 

0 

JWU 

- — -dy=  Cross  Product  Momentum  Thickness 

0 

"c  \N^ 

5,=  - — afy=  Crossflow  Momentum  Thickness 

7  J  112 


Cross  Product  Momentum  Thickness 


J  //2 
0  '■'e 
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Figure  1 .  Measurement  locations  and  definitions  of  i 


2: 


o 

di 

o 


-0.6  -0.5  -0.4  -0.3  -0.2  -0.1 


0.0 


Level  u 

F  1.43E0 

E  1.33E0 

D  1.24E0 

C  1.14E0 

B  1.05E0 

A  9.51  E-1 

9  8.56E-1 

8  7.61  E-1 

7  6.65E-1 

6  5.70E-1 

5  4.75E-1 

4  3.80E-1 

3  2.85E-1 

2  1.90E-1 

1  9.51  E-2 


s/t 


Figure  2.  U/U^g^  velocity  magnitude  contours  and  secondary  flow  streamlines. 


V 

3.41  E-2 

1.17E-2 

-1.07E-2 

-3.30E-2 

-5.54E-2 

.7.78E-2 

•1.00E-1 

-1.22E-1 

-1.45E-1 

-1.67E-1 

-1.90E-1 

-2.12E-1 

-2.34E-1 

-2.57E-1 

-2.79E-1 


Figure  3.  V/U,g,  velocity  magnitude  contours  and  secondary  flow  streamlines. 
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iog„(y/i)  S  '“S.ofy'') 


2.16E-1 
1.73E-1 
1.29E-1 
8.53E-2 
4.17E-2 
-1.96E-3 
-4.66E-2 
-8.93E-2 
-1.33E-1 
-1.77E-1 
-2.20E-1 
-2.64E-1 
-3.07E-1 
-3.51  E-1 
-3.95E-1 


velocity  magnitude  contours  and  secondary  flow  streamlines. 


Level  ii2 
F  3.51  E-2 

E  327E-2 

D  3.04E-2 

C  2.80E-2 

B  2.57E-2 

A  2.34E-2 

9  2.10E-2 

8  1.87E-2 

7  1.64E-2 

6  1.40E-2 

5  1.17E-2 

4  9.35E-3 

3  7.01  E-3 

2  4.67E-3 

1  2.34E-3 


sA 


normal  stress  contours  and  secondary  flow  streamlines. 
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iog,.(y«)  s  i°9.o(y'') 


Level  v2 


3.22E-2 
3.01  E-2 
2.79E-2 
2.58E-2 
2.36E-2 
2.15E-2 
1.93E-2 
1.72E-2 
1.50E*2 
1.29E-2 
1.07E-2 
8.60E<3 
6.45E-3 
4.30E-3 
2.15E-3 


ure  6. 


normal  stress  contours  and  secondary  flow  streamlines. 
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normal  stress  contours  and  secondary  flow  streamlines. 
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Figure  8.  i/WU,,,*  shear  stress  contours  and  secondary  flow  streamlines. 
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Figure  9.  shear  stress  contours  and  secondary  flow  streamlines. 
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Figure  11b.  Probability  density  function  of  v  fluctuating  velocity  at  s/t=-0.247. 
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Figure  11c.  Probability  density  function  of  w  fluctuating  velocity  at  sA=-0.247. 
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jre  12a.  Probability  density  function  of  u  fluctuating  velocity  at  sA=-0.291. 
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Figure  12b.  Probability  density  function  of  v  fluctuating  velocity  at  s/t=:-0.291. 
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Figure  12c.  Probability  density  function  of  w  fluctuating  velocity  at  s/t=-0.291. 
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Figure  17.  triple  product  contours  and  secondary  flow  streamlines. 
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Figure  20.  triple  product  contours  and  secondary  flow  streamlines. 
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Figure  21.  triple  product  contours  and  secondary  flow  streamlines. 
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Figure  22.  triple  product  contours  and  secondary  flow  streamlines. 
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Figure  23.  viI^/U,g®  triple  product  contours  and  secondary  flow  streamlines. 
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Figure  24.  uvwl\Jj  triple  product  contours  and  secondary  flow  streamlines. 
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Figure  25.  Flow  angle  (degrees)  contours  and  secondary  flow  streamlines. 
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jre  26.  Flow  gradient  angle  (degrees)  contours  and  secondary  flow  streamlines. 
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Figure  27.  Shear  stress  angle  (degrees)  contours  and  secondary  flow  streamlines. 
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Figure  28. 
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vorticity  vector  component  contours  and  secondary  flow  streamlines. 


Figure  29.  Turbulent  kinetic  energy/Ufe*  contours  and  secondary  flow  streamlines. 
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Figure  30.  transport  velocity  of  turbulent  kinetic  energy  contours  and  secondary  flow  streamlines. 
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Figure  31.  Sg  =1/^/7/^  structural  parameter  contours  and  secondary  flow  streamlines. 
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Figure  32.  Skewness  of  u  structural  parameter  contours  and  secondary  flow  streamlines. 
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•  Figure  34.  Skewness  of  w  structural  parameter  contours  and  secondary  flow  streamlines. 
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Figure  35.  Shear  stress  magnitude/Ure,^  contours  and  secondary  flow  streamlines. 
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jre  36.  Townsend’s  structural  parameter  contours  and  secondary  flow  streamlines. 
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Structural  parameter  contours  and  secondary  flow  streamlines. 
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#  Figure  38.  Rotta’s  T  eddy  viscosity  ratio. 
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Appendix  A*  Line  Plots  of  Flow 
Variables  and  Some  Structural 
Parameters 
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In  this  Appendix  line  plots  of  the  flow  variables  and  some  structural  parameters  are  given 
in  the  22.4®  coordinate  system.  Mean  velocities  (pages  53-60)  and  Reynolds  stresses  (pages  SI¬ 
TS)  are  plotted  separately  for  each  station.  For  triple  products  and  derived  quantities  the  stations 
are  marked  on  the  figures.  The  order  of  the  stations  is:  station  5,  a,  b,  c,  d,  e,  tb,  tc,  r,  s,  t,  u,  v, 
w,  X.  The  derived  quantity  and  structural  parameter  figures  are  presented  on  pages  87-104  after 
the  triple  product  figures  (pages  77-86). 

The  derived  quantities  and  the  structural  parameters  presented  are: 

Flow  angle  (FA),  page  87; 

Flow  gradient  angle  (FGA),  page  88; 

Shear  stress  angle  (SSA),  page  89; 

Mean  velocity  vector  magnitude  (Q),  page  90; 

Mixing  length  (L,„),  page  91; 

Turbulent  kinetic  energy  (TKE),  page  92; 

Townsend’s  structural  parameter  (A,),  page  93; 

v^/|t/p|  ,  page  94; 

Uv/TF,  page  95; 

(lP+M^)/v^,  pages  96-97; 

The  ratio  of  viv  stress  production  to  that  of  Uw  stress  production  C=  ^  tan(FGA)  gg. 

uv\an{FGA)  +  VW 

u  skewness,  page  99; 

V  skewness,  page  100; 
w  skewness,  page  101; 

Turbulent  kinetic  energy  transport  velocity  page  102; 

Longitudinal  axis  vorticity  component,  page  103; 

Shear  stress  magnitude,  page  104. 
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